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Linear superpositions of macroscopically distinct quantum states (sometimes also called 
Schrodinger cat states) are usually almost immediately reduced to a statistical mixture if exposed to 
the dephasing influence of a dissipative environment. Couplings to the environment with a certain 
symmetry can lead to slow decoherence, however. We give specific examples of slowly decoher- 
ing Schrodinger cat states in a realistic quantum optical system and discuss how they might be 
constructed experimentally. 

'd\- 

G^ ' I. INTRODUCTION 

a^ . 

^ : A Schrodinger cat state is a linear superposition of two quantum states that differ on a macroscopic scale. Whereas 
superpositions of quantum states are commonplace in the microscopic world, a superposition of macroscopically dis- 
tinct states is practically never observed. The understandingof why this is so has evolved considerably with the 

• development of the quantum mechanics of dissipative systems . Dissipative systems are systems that are coupled 
^-H \ to an environment (also called "heat bath" ) with a large number of degrees of freedom. The coupling allows for an 

. exchange of energy between system and bath. Therefore the motion of the system will slow down till only thermal 
fluctuations in equilibrium with the heat bath are left. This effect is also present in classical mechanics and takes place 
" . on a classical time scale, Tdass- But the coupling to the environment also very rapidly destroys quantum mechanical 
^ I ' interference effects, and leads to an effectively classical behavior. The decoherence arises on a time scale Tdec which 

, is typically much shorter than Tdass- As a general rule, if a seperation D in phase space can be assigned to the states 

• involved in the Schrodinger cat state, the two time scales behave as TdecjTciass ~ i^/DY where r is some positive 
CO ' power depending on the details of the coupling and the environment P-p[. This separation of time scales has also 

• been termed "accelerated decoherence" . The validity of this picture starts to be confirmed in experiments of Haroche 
et at, who were able to turn on the coupling to an environment in a controlled way and measure the decay of the 

^ \ coherences §,|l0|. 

Qh. The basis in which the density matrix becomes diagonal depends on the coupling to the environment, as was pointed 
4!) ' out by Zurek Suppose that the coupling Hamiltonian contains a "coupling agent" A, Hint = ^/({x, p}), where 
^ _ /({x, p}) is some function of the environmental coordinates. Accelerated decoherence arises in the basis formed by 
^ ' the eigenstates of A, the socalled pointer basis. Thus, after the short time Tdec the reduced density matrix of the 
system will be diagonal in this basis, but may still contain off-diagonal matrix elements (coherences) in another. 
' It has recently been shown, however, that even in the pointer basis accelerated decoherence can be absent if the cou- 
.5^ I pling to the environment has a certain symmetry ||^,[ll|-|l5|| . Indeed, suppose that AjV'i) = a\ipi) and A\iIj2) — 01-02), 
■ i.e. a is a degenerate eigenvalue of A, then the environment cannot distinguish between states jV'i) and |'(/'2)- A 
Jh [ Schrodinger cat state = ci\ip) + C2IV') is therefore stable against dephasing. It may loose its coherence on the time 
. - - ■ scale on which the symmetrical arrangement of the state is destroyed, which can be as large as Tdass- 

In this paper we give specific examples of such longlived Schrodinger cat states in the context of superradiance with 
a well known, experimentally verfied damping mechanism. We calculate the decoherence rates and propose a way 
how the Schrodinger cat states might be realized experimentally. 



II. SUPERRADIANCE AND THE SCHRODINGER CAT STATES 

The system that we consider is well known from the phenomenon of superradiance : A large number of identical 
two-level atoms in a cavity radiate collectively. We assign formally to each atom a spin-i operator S^*-* whose z- 
components tells us whether the atom is in its lower or upper state. If all atoms couple in the same way to a resonant 
electromagnetic field mode in the cavity, the system as a whole can be described by a collective variable J — X^t^i S*-*-* , 
the so called Bloch vector [N is the number of atoms). It can be thought of formally as an angular momentum. The 
coupling of the atoms to the resonant electromagnetic field mode (with creation and annihilation operators and a) 
be of the form 
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Hint = TT-g{J+a + J-o^) , (1) 

with some coupling constant g. Dissipation ultimately arises by the cavity's non-ideal mirrors, i.e. the field mode 
is itself damped and photons may leak out of the cavity at a rate k. Under the assumption of weak coupling 
{g^/N / n <^ 1), low temperatures {ksT <C Tiui, the level seperation of the atoms), and the Markov assumption 
(t ^ I/'*) the master equation 

^Pir) = Y^{[J-,p{r)J+] + [J-p(r), J+]) (2) 

for the reduced density operator p was derived in . We have written the time t dimensionless in units of the inverse 
classical damping rate as r = ^ig^tjK^ i.e. t is already in units of Tdass- The classical limit is obtained by j ^ oo 
and yields a spherical phase space in which the orientation of the classical Bloch vector can be parametrized by the 
polar angles 9 (reckoned against the Jz-axis) and ip (the azimuth of the equatorial projection reckoned against the 
Ja;-axis). The variables p — Jz/j = cos 9 and (p aquire the role of classical canonical momentum and coordinate [ p^ . 
Since the classical phase space contains 2j + 1 states, we may think of Ti as scaling like 1/j. The Bloch vector behaves 
classically like an overdamped pendulum, d9/dT = sin 9 and (p = const.. Many of the consequences of (|^) have been 
confirmed experimentally |19|| . 

The states that correspond as much as possible to classical states are the so called angular momentum coherent 
states I7) = \9,ip) [^|T]. They correspond to a classical angular momentum pointing in the direction given by 
the polar angles 9 and ip, with the complex label 7 given by 7 = taii(9/2)e^'^ . They have minimal uncertainty, 
ApAq ~ In terms of \ jm) states one has the expansion 




|7> = (1 + 11*)-' E ^'""l/ ■ _ 1^'"^) • (3) 



Coherent states may be more familiar from the harmonic oscillator, where they are eigenstates of the annihilation 
operator. The compactness of Hilbert space prevents the existence of exact eigenstates of J_, but one can show that 
the angular momentum coherent states are approximate eigenstates of J_ in the sense that the angle between J-\j) 
and 1 7) is of the order of They therefore qualify as pointer states in the limit j 00. 
In the following we will consider the fate of a Schrodinger cat state 

W^^i\ll)+\l2)), (4) 

where Af is the appropriate normalization constant. Note that the two components can indeed be macroscopically 
distinct if the number of atoms in the superradiance experiment is large. We will show that the damping (H) leads 
in general to accelerated decoherence. Our central prediction is, however, that accelerated decoherence is absent for 
Schrodinger cat states with 7172 = 1- Such two states correspond to two classical angular momenta arranged sym- 
metrically with respect to the equator 9 — 7r/2 (i.e. 92 — it — 9i) on a great circle (^i = </'2 = const.. They decohere 
only on the classical time scale, T^ec ~ Tdass- 

The initial density matrix corresponding to the state (||) reads /5(0) = |A/'P(|7i) (71I + |7i)(72| + |72)(7i| + |72)(72|)- 
Since the evolution equation of /5(t) is linear it suffices to concentrate on one off-diagonal part /5(r) evolving from 
p(0) = |7i)(72|- Its decay is conveniently studied in terms of one of the norms 

Ni{t) ^trpp^OT (5) 



A^2(t)= V \Pm,m,\, (6) 



■mi ,m2 



where Pmim2 — (j"T'i|pIj"^2)- Both norms are obviously zero if coherences in the 7 basis are absent. The time on 
which they decay to zero sets the decoherence time scale Tdec 

Our prediction of slow decoherence, Tdec ^ Tdass for 7172 = 1, is based on three analytical results: 



1. The initial time derivative of Ni{t) reads 

1= — 2j (sin"^ 6*1 + sin"' 6*2 — 2 cos((y92 — ((Ci) sin6'i sin6'2) (7) 



dT '^=" 

- ((l + cos6'i)2(l + cos6i2)^) 
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Thus, for ip2 = </?i and sin 9i = sin 02, the term proportional to j vanishes. This marks the absence of accelerated 
decoherence. 

In the following we will restrict ourselves to ipi — Lp2 — Q and denote 7.; = |7i| throughout the rest of the paper. 

2. In the particular case of the states I71) — and I72) — \i — j) (corresponding to 71 = and 72 = 00), the 
exact time evolution N2{t) is given by N2{t) = e""^ for all times! The coherence decays on the classical time 
scale, even though the two states are macroscopically as distinct as possible. 

3. For finite times r with jt ^ 1, a semiclassical evaluation of the norm N2{t) for Lpi = if2 — Q leads to 

N2{t) ( ^. (71-72)^(1-7172)^ ^ , ^o,-,/.^^ 

^ ^ ((i+7f)(i+7i))^ V + ""^"'^^ ■ 

This means accelerated decoherence as long as 71 7^ 72 and 7172 7^ 1. If, however, 7172 = 1 then the next order 
in 1/j shows that 

iV2(0)-°"Pl^ U + lj 2(7? + 1)4 " ) 

•(l + 0(l/j)). (10) 

The expression in the exponent is correct up to and including order (jV)^. Obviously, accelerated decoherence 
is absent for 7172 = 1. Indeed, a single coherent state 71 = 72 = 7 leads, in linear order, to almost the same 
decay, 

'^^^")-expl-7M^) -1 ■ (11) 




iV2(0) 

It is clear that the symmetry of the sine function under 9 ^ w — 9 underlies the longevity. As mentioned the 
coherent states I7) are approximate eigenstates of the coupling agent J_ and therefore qualify as pointer states. The 
symmetry of the sine function means that the corresponding "approximate eigenvalue", jsin^e^*'^, is doubly degen- 
erate. A linear combination of two vectors of the pertaining subspace is stable against decoherence. We have here 
the interesting situation that the deviation of order 1 /j from exact degeneracy is small enough to compensate for the 
acceleration factor j of the decoherence rates. Similiar but more general statements about the decoherence in the 
presence of degeneracy breaking coupling agents can be found in |1 1| . 

We now describe in some detail the derivation of the three mentioned results. The first, eq.(H) can be obtained 

straight forwardly by inserting the equation of motion for p(t) into 7^^^—^ = tr(^p^ + p^^)- 

The second result is an exact solution of the master equation (|^) ]l7|p^ . To see this we write (||) in the |jm)-basis 
and introduce mean and relative quantum numbers, m — "1+™^ and k — ™^~™^ . The latter quantum number is 
conserved and therefore only enters as a parameter, Pmi,m2 — Pm{k,T). The solution can be written as 

3-\k\ 

Pmik,T)^ J2 £'m„(fc,T)p„(fc,0) (12) 

n=-j + \k\ 

in terms of the dissipative propagator Dmn.{k,T). With the help of the abbreviations Qmn — \j+Tn)l(j^)\ ^^'^ 
5i = j(i + 1) ^ — 1) 0116 exact Laplace integral representation for Dmn{k, r), 

"^Z Q'rn — k,n — kQm+k,n+k ^ ^ 



27ri 



/ d.e"/^n— ^2' (13) 



where the real parameter b should be chosen larger than the largest pole in the denominator. Depending on m and 
n, in general a large number of poles contribute to the integral. But for the case at our interest (mi = j, m2 = — j), 
we have m = = n, fc = j and therefore only one pole contributes. We immediately obtain N2{t) = pa{j,T) = 
Doo{j,T)po{j,0) = exp(-r). 

Starting point for the third result is the short time propagator 



3 



(n-TOj! \j J 

derived in It is valid up to times (ji _ ttj),- <^ Since we are going to sum over m and n, both m + n 

and m — n can be of order j and one should therefore have jr <C 1. The expansion (^) shows that (at ip = Q) all 
initial density matrix elements are real and positiv. Since also the short time propagator is real and positiv, the same 
is true for the density matrix elements to the later time r, such that the norm A'2 simplifies to 

^^2(r) =^p™(fc,r) = ^ i?„„(fc,T)p„(fc,0). (15) 

fe,m m,n,k 

Note that in general iV2(0) 7^ 1, such that in order to study the decay of the coherences one should look at the quantity 
n(T) = 7V2(r)/7V2(0). 

The first step in calculating n(r) is to evaluate the sum over the final states m. This summation does not envolve the 
initial density matrix at all, so we can define S{n,k,T) = J2m ^mnik^r) and N{t) — ]^S{n,k,T)pn{kTO). The 
latter two summations will eventually be done by Laplace's method, that is by an asymptotic expansion in l/j for 
j — > 00. In order to remain consistent with the limitation jr <C 1 we will have to rescale time as 

"J («, 

with f <C 1 kept fixed for j — > 00. This is indeed reasonable since we expect decoherence of an ordinary Schrodinger 
cat state on a time scale l/j, that is N{t) ~ -/V(O) exp(— jr). If we kept t fixed for j ^ 00 we could only recover 
N(t) = 0, whereas we can study the decay if we keep f fixed. The precise calculation of S{n, fc, t) by summing over 
m is difficult. However, since we shall be interested only in short times (jr <gC 1) we can reconstruct >S'(n, k, r) from 
its initial time derivatives. Let us write the derivatives directly in reduced coordinates that become continous in the 
limit j — + CXI, v = n/j and 77 ~ k/ j, as 

9^(n,fc,r=f) 1 



d^Sin,k,T=f) _ f dS{n,k,T=j) \ 



Wi = 



2 



+ Wi(w2 — Wl) (17) 



f=0 
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For the further derivation we will keep S in the form of a polynomial, S{iy, rj, f) = l+a{v, rj, \/ j)f+h{v, rj, 1/ j)f'^-\-0{T^) 
with a{h',r], l/j) = dS/df and b{u,r], l/j) = ^d^S/dr^ \f=o- Higher terms in f can be included, but already the 
first two terms account rather well for the initial decay of the coherences. 



We now transform the remaining two sums in N2{t) (eq.([T^)) into integrals by Euler-Maclaurin summation and 
finally integrate via Laplace's method. To this end we first write the coefficients /0„(fc, 0) in the expansion of the initial 
density matrix as continuous functions of the reduced coordinates r/ and v = n/j. The result can be cast in the form 



Pr^ik.Q) = Ce^^-^''^^\ (18) 



where the "action" 6*0 is given by 



So{v, ry) - (1 - u) ln(7i72) + 77 In ^ - Up{l - - v) (19) 

71 2 

+p{l + i' + T])+p{l-v + 7j)+p{l + iy- T])) (20) 

p{x) = xlnx . (21) 
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The action is correct up to lowest order Terms of order 1/j have been absorbed in the prefactor C which turns 
out to be independent of v^rj and therefore wiU cancel in ^(t). We thus have to order 



N2{T)^2fC J dT^J diy {1 + ai,,,7j,l/j)f + bi,,,r^,l/j)T^) eMjSo{'^,v)) ■ (22) 

The prefactor 2 stems from the fact that k and n can be simultaneously integer or semi-integer, but cancels of course 
as well in n(r). Since the slow decoherence is a 1/j effect, a carefull asymptotic expansion of the double integral 
in powers of 1/j is in order. First of all, we expand the coefficients a and b as a{v,r],l/j) = Y^'^=o'^k{'^,v)jk and 
77, 1/j) = X^fc^o ^^j^' "^^^ ^^^^ terms in the expansion read {w = ^{1 — {1/ — rj)^) {1 — {v + r/)^)) 

ao = w + 7?^ - 1 + i/^ (23) 
l-iy^+ij'^-w 

ai = V (24) 

w 

bo = ^al (25) 
bi = (aofli - + w-r]'^ - 1)) (26) 
62 = 77^ f (2«oa2 + af){{rj^ - u'^f - 2{Tf + u'') + l) - 2 

2Tf - IQv^rf + 8j/^77* - IQv^ - Si^^rj^ + 2r]^ + 8j/^ + 2tj^^ 
+^ (Sr]^ + 7iy^ - IQv'^rf - lOi/^ - Qif + 3) . (27) 

For the expansion to make sense w should be sufficiently far (^ 1/j) away from 0. The coefficient 02 will not be 
needed. If we insert the expansions in the integral (E2[) we encounter functionals of the type 



2w2 



/ dr^ / dyf{y,7i)e^^-^''^^\ (28) 



where f{v,Ti) can be 1 or any of the coefficients in ( ^3|) to (^7|). The large parameter j in the exponent suggests to 
integrate by the two dimensional Laplace method. To obtain correctly the first few orders in 1/j for n(r) one needs 
an extension of Laplace's method to higher orders, i.e. we also have to expand /[/] in powers of 1/j'. 

m = [/] + [/] + ^/^'^ [/] + .... (29) 

The higher orders can be obtained from systematically reexpanding in 1/j the following exact representation of /[/] 
11, 

/[/] =e^-^°('^-*)^|dct(^)ri^-— ^(L5)' (/(^,r?)e^-«(''''')) , (30) 

where (i/q, 770) denotes the position of the maximum of Sq. The matrix cr contains the four second derivatives of 5 at 
the maximum, cr^^jj — d^dr/So the new "action" i?(i^, 77) is the deviation of Sq from its quadratic approximation, 

Rii., rj) EE 5o(^, V) - ^0(^0, %) - ^ ' (31) 



2 XV^VoJ 

and Ls is the second order homogeneous differential operator 



Ls = {-<r-^(^:], (':))■ (32) 



dr,J ' \dr, 



The angular brackets denote a scalar product. It is assumed that the integration range contains exactly one maximum 
with d^So{i', rj) — = (9,,S'o(i^, 77). Due to the construction of R, all derivatives of R of lower than third order vanish 
and the term Lg{fe^^) in (|3^) is therefore a polynomial in 1/j at most of the order 2//3. In order to get the first 
order in 1/j beyond the usual Laplace method (the term with / = 0) one has therefore to go up to Z = 3. 
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Combining the expansions of a and b in terms of ak,bk and the expansions of the corresponding functionals I[ak] 
and I[bk] according to (30) we find 



n(T) — I + f 



+/W[ao])/(°)[l])+o(l 
I^°Hbo] 1 



(/("'[ao]/(i)[l]-(/(°)[ai] 



/(o)[i] ,(/(°)[i])^"^''°'f'°^''''™ " (/(°)[6,] + /(^)N)/(°)[i]) 



j2(/(0)[l]) 

[!]/(") [fei] + /(I) [60])/^''^ [1] - [60] (/('^ [1]) 



The point of maximum So is located in our case at 



m = 



ll - ll 



1 - ihl 



(l + 7?)(l + 7l) 



(l + 7?)(l+7l)' 



(33) 



(34) 



To proceed further we distinguish two different cases. 
Case 1: 71 7^ 72 and 7172 ^ 1 

Here aheady the ordinary Laplace method leads to a meaningful result, 

/(°)[ao] f , 2(71 -72)^(1 -7172)' 



/(o)[l] 



((l+7?)(l+7l))' 



/(o)[i] ^ M'^o,Vq) = 2^0*^^°''^°^ 



(35) 
(36) 



We have therefore ^(t) = 1 + 00(1^0, '?o)7' + 500(^01 '7o)^ j and thus, correct up to 0{f^) the anounced result (g), where 
we have resubstituted f in terms of jr. 

Case 2: 7172 = 1 

The leading terms (^5|), (|3^ ) now vanish due to ao(t'o, ?yo) = 0. From the prefactor of the term linear in f in ( ^3|) only 



i V7f + 1 



survives, from the quadratic term only 



l /W[ao] _ 1^7f-l 
1 /("'[^'2] 



j2 /(0)[i] 4f. 



(37) 



(38) 



as the reader might verify in a straightforward but lengthy calculation. In particular, the coefficient proportional to 
1/j in the term is zero, such that to quadratic order in f ^(t) depends only on r, not on jr: 



n(T) = 1 



(1! 


1— ' 






+ 1. 





9 2 

72 - 7i 



(l+7?)(l+7l) 



+ 1 T2+0(f3) 



(39) 



The agreement of Eq. ( |39[ ) with exact numerical results extends well beyond the range r ^ 1/j for which the theory was 
made initially. This is not surprising since the main contribution to the norm comes from a region where n ~ ~ to, 
so that the limitation on the validity of the short time propagator is much less severe than in the asumed worst case 
where both n — m and n + m are of order j. The agreement with the numerical result becomes even better and leads 
to rather precise results even for t ~ 1 if we rewrite the decay again in exponential form as was done in (0). 
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III. POSSIBLE EXPERIMENTAL REALIZATION 



It has been experimentall y v erified that (|) describes adequately the radiation by identical atoms resonantly coupled 
to a leaky resonator mode p9[ in a suitable parameter regime (see the discussion in the context of the equation). It 
should therefore be possible to observe the slow decoherence of the special Schrodinger cat states. We now propose 
a scheme for their preparation. Starting by all atoms in the ground state and with the field mode in its vacuum 
state, a resonant laser pulse brings the Bloch vector to a coherent state \9,Lp). Note that the cavity may be strongly 
detuned (with respect to the atomic transition frequency) during the whole preparation process (detuning 6 ^ k). 
The dissipation mechanism (H) is hereby turned off and the system evolves unitarily with a Hamiltonian containing 
a non-linear term oc {g^ / 5)J+J-. The free evolution during a suitable time will split the coherent state |0, i^) in a 
superposition of \9, ip') and 16*, ip' + tt) as described in Finally, a resonant n/2 pulse brings the superposition to 
the desired orientation symmetric to the equator, by rotation through the angle 7r/2 about an axis perpendicular to 
the plane defined by the directions of the two coherent states produced by the free evolution. At this point the cavity 
can be tuned to resonance, thereby switching on the dissipation mechanism and one can study the decay of coherence. 



IV. CONLUSION 



We have shown that a certain symmetry of the coupling to the environment leads in the phenomen of superradiance 
to the existence of longlived coherences of superpositions of macroscopically distinct quantum states. Even though the 
components of the linear superposition are not exact eigenstates of the coupling operator to a degenerate eigenvalue, 
the deviation from degeneracy is small enough for the coherences to decay on a classical time scale only. We have 
proposed a preparation scheme with which such Schrodinger cat states might be realizable experimentally. 
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